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^ : Abstract 

J2 ■ We consider the adsorption of a random heteropolymer onto an interface 

within the model by Garel et al. |jl] by taking into account self-interactions 
, between the monomers. Within the replica trick and by using a self-consistent 

^ I preaveraging procedure we map the adsorption problem onto the problem of 

] binding state of a quantum mechanical Hamiltonian. The analysis of the 

latter is treated within the variational method based on the 2-nd Legendre 
transform. We have found that self-interactions favor the localization. The 
, effect is intensified with decrease of the temperature. Within a model without 

lO i taking into account the repulsive ternary monomer-monomer interactions we 

predict a reentrant localization transition for large values of the asymmetry 
of the heteropolymer and at low enough temperatures. 

O '. PACS numbers: 61.41. -he, 05A0.+j, 03.65.-w 
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INTRODUCTION 



The behavior of heteropolymers at interface between two immiscible (incompatible) sol- 
vents has been intensively studied in recent years |l|]- [|16| since it has an obvious importance 
in biological applications (proteins and membranes) and application in different fields of 
industry such as biosensor, pattern recognition applications, glues, paints etc. [|18[. Experi- 



ments |19[, |22] and numerical simulations I^Tl, |23[ have shown that a random heteropolymer 



may localize, reinforcing the interface between two incompatible solvents and reduce inter- 
facial tension. Recent theoretical efforts have been devoted to understand the fundamental 
physical mechanism governing the localization of a random copolymer onto an interface 

In the simple model introduced by Garel et al. |I| only the interaction of the monomers 
with the solvent, but not the self- interactions between the monomers were taken into ac- 
count. Considering a. A — B copolymer at the A — B interface A monomers prefer to be 
in the A-half-plane while B monomers prefer to be in i?-half-plane. Obviously there is a 
sort of frustration in such system because the complete separation of A monomers into one- 
half and B monomers into another half plane is forbidden by polymer bonds. The analysis 
performed in |]l| showed that the localized state of a random heteropolymer chain in the 
presence of a selective interface can be imagined as consisting of blobs with majority of A 
monomers and minority of B monomers in the A-half-plane and vice versa for the blobs in 
the 5-half-plane. It was shown that random copolymer always localizes for statistically sym- 
metric heteropolymer, whereas a delocalization transition was found if the heteropolymer is 
asymmetric. The heterogeneity in the chemical structure of the polymer, which results in 
self-interactions between the monomers, may have a considerable impact on their bulk ther- 
modynamic behavior which consists in a segregation into A- rich and 5-rich domains. 
In the case of a single heteropolymer in solvent the self-interactions favor the collapse of 
the heteropolymer ||2^. The self-interactions between monomers may play an important 
role for adsorption as it was noticed in [Q, [^. However, the influence of self- interactions 
on the localization behavior of the random heteropolymer is poorly understood. The effect 
of the excluded volume in the adsorption of a heteropolymer was recently investigated in 
| 26| . In the present paper we will consider the influence of the direct monomer- monomer 



interactions on the localization behavior of a random heteropolymer onto an interface. 

The article is organized as follows. In Section II we will introduce to the model and the 
formalism we use. Section III contains our results and Section IV conclusions. 



I. MODEL AND FORMALISM 

We study the problem of a random heteropolymer by using the model proposed in [|l|. 
Let us consider a long two-letter A — B heteropolymer chain in the presence of the A' — B' 
interface between two incompatible solvents. Alternation of different type of monomers 
along the chain is assumed to be random, and each monomer is assumed to interact with 
an external potential, which takes positive and negative values depending on the position 
with respect to the interface. The partition function of the polymer in the presence of the 
interface is given by the following Edwards functional integral over the trajectories r(s) 
(0 < s < A^) of the polymer chain 
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^{C(s)} = y"/^r(s)exp|-^^ ds(^^^ + w ds as)sgn[z{s)] + 

ds ds' as)as)K<s)-v{s'))^, (1) 

where a is the Kuhnian segment length, z{s) is the Cartesian component of r(s) transversal 
to the interface, w and xo are measured in units of ksT. The first term in the exponential 
of (|I]) corresponds to the elastic energy of the polymer chain, the second one describes the 
monomer interaction with the medium, which is governed by the random parameter C,{s), 
which is assumed to be Gaussian distributed with the distribution function 



P{C(s)} cx exp 



1 



2A, 



N 

{Cis)-Co)'ds 



Jo 



(2) 



The distribution function @ of the random "charges" ({s) is completely characterized by 
its two moments, {({s)) = (o with (q related to the asymmetry in the composition of the 
copolymer, and {C{s)C{s')) = (q+Ao6{s — s') where Aq being the variance of the distribution. 
The interaction potential is chosen as a step function. Such a choice is legitimate, if the 
interface width is much smaller than the size of an average excursion loop. The last term 
in Eq.(|l]) describes the short-rang interaction between monomers, where Xo is the effective 
interaction potential (the second virial coefficient). The sign of the last term in (|I]) is chosen 
so that for xo > the like monomers attract each other while the unlike monomers repel each 
other. The model described by Eq.(0) admit that the self-interactions occur independent 
of whether the monomers are on the left or on the right side of the interface. This may be 
the case if even the favorable solvent for a given type of monomer is slightly poor. If the 
favorable solvent is comprised of the same monomers, then the self-interactions between the 
like monomers in their own medium are expected to be zero. The generalization of the model 
where the self-interactions are switched off, if they are in their own medium is possible. The 
effect of ternary interactions, which will prevent the collapse of the blobs of the adsorbed 
heteropolymer due to attractive self-interactions, will be considered in a separate paper. 

The random "charges" ({s) in (|l]) are considered as quenched variables, so that in or- 
der to obtain the free energy one has to average ln(Z) over all possible realizations of 
monomer sequences. For this purpose we use the replica trick consisting in introduction 
of n copies of the system with the same quenched variables C{s), and using the identity 
ln(Z) = \imn^o{Z^ — l)/n in averaging over C{s). Thus, at the intermediate stage we con- 
sider the average (Z^) where (...) means average with the distribution function (Q). The 
partition function can be written as 



/n 
Y[Dr^{s)exp{-Ho-Hint}, 
a=l 



(3) 



where 
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Hint = -W 



/ rfsC(s) Vsgn[r,Gs)]-Xo / ds ds' Cis)as')y25 Ms) - r^is')) . 

(5) 

To average over ({s) in (|^) we expand in Taylor series in powers of Hint, carry out the 
average, and reexponentiate the obtained expression. This results in writing the result of 
the average in the exponential as 

In ((exp {-H.r^t})) = - {H.nt) + \ {{HL) - {H^nt? + ...)• (6) 



Restricting the expansion in @ to quadratic terms in Hint we obtain 



(^") = y"ni5r,(s)exp|-^ Cnds^ 



(7) 



with 

^ a=l ^ a=l a=l 

-w'^Aq ^ sgn(r„(s))sgn(r^(s)) - 2u;xoCoAo ^ sgn(r„(s)) / 6 {r,3{s) - r,3{s')) ds' 

a, 13=1 a, 13=1 "^^ 

" pN nN 

2x^Co'AoE/ / ds"5{v^{s)-r^{s'))5{v,{s)-v,{s")) - 

n „jv 

XoAo E / ^ (r-(^) - ^ (r/3(.) - r,(.')) ds' . 

Notice that contains more than one integration over the contour length. Due to this it 
is not possible to reduce to a quantum mechanical Hamiltonian as it is the case, if only 
a single integration over the contour length appears (see for example 0). In the following 
we will preaverage (^) in order to reduce it to a quantum mechanical Hamiltonian. For this 
purpose let us consider the quantities 



p(r,)= / 5(r,-r,(s')) ds' , (9) 
Jo 

pN 

g(r^, T/s) = 6 {to, - r^is')) 6 (r^ - r^(s')) ds', (10) 
Jo 

where p{ra) is the microscopic density associated with the polymer and g(rQ,,r^) is the 
replica overlap parameter [^. We will preaverage p{Ta) and g(rQ,r^) in Eq.(^ over ra{s') 
and rfs{s') according to 

(p(r„)) = ds j dZ J dZ'G{Z,N;z,s)G{z,s;Z',0) x 

jdR^J di?:,Go(i?,„ N- r,„ s)Go(r,„ s; R'„, 0) (11) 



and similar for q{ra,rf3). Go{R„, N;r,„ s) in ([TlD is the free Green's function of the longi- 
tudinal degrees of freedom of the polymer. Notice that the average over the longitudinal 
coordinates of the polymer with the unperturbed Green's function is an approximation. We 
use the spectral representation of the transversal Green's function in Eq.([Tl|) G{z, N; z' ,0) 
through the eigenf unctions ipk{z) 

G{z, z'; N) = J2 e-'^'-Mmi^') ■ (12) 

k 

The Green's function G{z, z'\ N) satisfies for > the differential equation 

~dN-~j^'^^~r^' (13) 

which is remarkably similar to the Schrodinger equation for a quantum particle in an external 
potential 

ot 2m 

The mapping of the latter onto Eq.(|TBp occurs by using the following replacements: t —>■ iN, 
h^T, a^T 3/m. 

In the case of a discrete spectrum with the energy of the ground state being negative, 
the main contribution in the series (0) for large originates from the ground state (ground 
state dominance) The preaveraging of p{ra) and q{ra, r^) given by Eqs. (|9|- |T0|) according 
to (|lT]) in the approximation of the ground state dominance gives 

(p(r„)) ^ a \iJko{zo.)\\ (14) 
(g(r„,r/3)) ^ a {tpkoiza)]"^ l^koizp)^ , (15) 



where a = N/S is an average monomer density per surface area pOf. We expect that the 
ratio cr is finite for large N and 5*. 

The substitution of (p^lf), (p!5|) into the Hamiltonian (|]) gives the effective replica Hamil- 
tonian as 

n n n 

Hn = -D^Vl- wCo ^ Sgn{Za) - XoCo^ ^ {-^{Za) f - 
a=l a=l a=l 

^ n n 

-w^Aq ^ sgn{za)sgn{zp) - 2wxoCoAocr ^ sgn{za) \^pizfs)f - 

a,f3=l a,p=l 
n 

(2CoV + Ao)x^Aoa ^ \^{z^)\'\^P{zp)\\ (16) 

Q,/3=l 

where we have introduced the notation D = Due to the average over the longitudinal 

degrees of freedom of the polymer by using the unperturbed Green's function, the problem 
becomes one dimensional. 

The investigation of the adsorption of a random heteropolymer chain is equivalent to 
the study of the ground state of the Hamiltonian if„ given by Eq . ([161) . Without taking into 
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account the self-interactions, the case which is obtained from Eq.(|T6D by putting xo = 0, 
Stepanow et al. [^] apphed a novel variational principle for Green's function associated 
with the Hamiltonian Hn. The latter generalizes the well-known Rayleigh-Ritz method 
in Quantum Mechanics for nonstationary states. The variational principle for the Green's 
function can be outlined as follows |^ . The start point is the Dyson equation for the Green's 
function G 

= -G-^ + Go^ + (17) 

which is considered as a stationarity condition 6F{G)/6G = of a functional F{G), which 
is found in a straightforward way as 

F{G) = -tr In(G') + tr Gq^G + tr (18) 

where the bare Green's function is defined as Gq^ = uj + Hq, with Hq being the unperturbed 
part of the Hamiltonian, is the interaction part of the Hamiltonian, and uj is Laplace 
conjugate to the chain's length A^. The functional F{G) given by Eq.(|T^ is the particular 
case of the generating functional of the 2-nd Legendre transform in Quantum Field Theory 
and Statistical Physics Notice that without preaveraging of p{Ta) and q{Ta,r/^) 

given by Eqs. (p|-p!OD the problem under consideration does not reduce to a quantum mechan- 
ical problem and instead of (|18|) we have to use the extremum principle associated with the 
second Legendre transform |^2|]. Assuming the ground state dominance we choose the 



n-replica trial Green's function as 

n 

G{kik2\ z, z'\ i) = W exp{-ekt)i'{kik2; Za)^{kik2; z'^), (19) 



where the trial wave function is chosen as 



HhM; ^a) = \ (e-'^^"^(^a) + e^'-^-d{-z^)) . (20) 

y ki + k2 

Notice that the energy Sk = — -D/cg is negative, and is a function of k2- Computing the terms 
in (p!8| ) by using (|19]), ( pOD gives the extremum functional as 

Tpr, , ^ ^ f , N , n{Dkik2-ek) nwCo f k2 - ki 
F{kik2) = m{uj + nSk) + 



ijj -|- n£fc ijj + nejfc \kx\k2 

n AqW^ n{n — 1) Aqw"^ fk2 — ki^'^ 



2 uj + nek 2 to + nek V ^1 + ^2 

50 



nXoCoCr kik2 n^wxoCoAoa 2kik2 f k2 - ki 



uo + nek ki + k2 uj + nek ki + k2 \ki + k2 
{n + 3n2) (2CoV + Ao)XoAofx / 2kik2 



(21) 



12 u; + nek + ^2, 

The stationarity conditions 

dF/dki = dF/dk2 = 0, (22) 
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give in the limit n = the following equations 



3kl + 9klk2 + 9kikl + 3kl + GCih + ^2) + 6A(A;i - ^2) - 

3xC^cTk2{ki + ^2) - 2(2CV + A)x^Aakikl = 0, (23) 

3klk2 + 9klkl + 9kikl + 3kl + 3C{kl - kl) - 6Akik2 - 

3xC(rkik2{h + k2) - (2CV + A)x^A(ykikl = 0, (24) 

where we have introduced new quantities A = Aqw^, Q = (qw, x = Xo/w'^ and have put 
D = 1. Solution of this system of equations gives us the localization lengths .^1 = 1/fci and 
.^2 = 1/^2; which describe the localization of the random heteropolymer onto the interface. 



II. RESULTS 

The present problem was studied in (see also 0- |]T^) without taking into account 
the self-interactions between the monomers. It was found in |P that the localization- 
delocalization transition occurs at the temperature = where the parameter ki be- 
comes zero and localization length ^1 becomes infinite. The value of T^. coincides exactly 
with that found in ^ by using a different method. The condition ki = Q means that the 
heteropolymer delocalizes in the right hand half plane (z > 0). 

Unfortunately, the equations (^), (|2^) cannot be solved analytically, so that we analyze 
them numerically. First we have examined the influence of the monomer-monomer inter- 
action on the localization lengths. In Fig. 1,2 we present the dependence of ki and ^2 as 
a functions of the asymmetry parameter Co for different values of the parameter xo- Fig-1 
corresponds to xo = and recovers the results of with values k\ = VqK/9 and A;^ = 2k^ 
for statistically symmetric heteropolymer, (^0 = 0. At the critical value of the asymmetry 
Q = 2A/3, the parameter ki becomes zero, so that the heteropolymer delocalizes in the 
right hand half-plane. The increase of the parameter xo (Fig-2) leads to the increase of the 
critical value of the asymmetry parameter at which the delocalization transition takes 
place. This means that the self-interactions favor the localization of the heteropolymer. 
This can be explained qualitatively as follows. The self-interactions being attractive results 
in a decrease of the size of the blobs on both sides of the interface, which result in a de- 
crease of their repulsion from the interface. The self-interactions influence to great extent 
the larger blobs, i.e. the blobs which are on the right hand side of the interface, while as 
in the case of the excluded volume interactions the effect of interactions is proportional to 
xVN with being the number of monomers in the blob. At given xo > and large enough 
Co and T (see Fig. 2, curve 2) the parameter ki becomes again positive that means that in 
this range of parameters the polymer is localized. Thus, as a function of the asymmetry 
parameter Co < Co < Co with Co and Co depending on x, there is window where the polymer 
is delocalized. For Co < Co ^^e self-interactions favor the localization driving the critical 
value Co to larger values. It is surprising that for Co > Co polymer localizes again. We 
interpret this as follows. We expect that in the model under consideration it is reasonable 
to impose a restriction on the maximal value of the asymmetry parameter C™'^^5 which is of 
order of Ao [EH], so that Co < C™^^- Although at Co = Co^^^ the heteropolymer is on average 
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homopolymer, due to the difference of the variance Aq from zero the typical polymer is still 
heterogeneous. The self-interactions between the monomers in Hamiltonian (16), which are 
attractive between the monomers of the same kind and repulsive between the monomers of 
different kind, have the tendency to cause a microphase separation between the A and B 
monomers at large values of asymmetry parameter (q. Then, the polymer can be imagined 
as consecutive units of microphase separated blobs, each of them prefers to be on the left or 
on the right with respect to the interface. This is expected to favor the localization of the 
heteropolymer. We expect that taking into account the incompressibility of the blobs will 
disfavor the reentrant transition. However, the prediction that the self-organization of the 
random heteropolymer due to self-interactions will favor the localization will remain valid. 

Notice that both ki and k2 increase with increase of Co and become equal at some (o 
(curve 1 in Fig. 2). For larger values of Co parameter ki becomes larger than k2. This signals 
that the blobs on the right hand side of the interface collapses to a size smaller than the blobs 
on the left hand side, which is reasonable due to the fact that the effect of self-interactions 
is amplified by the number of monomers in the blob. However, this region is beyond the 
physical realization of the present model, since we do not take into account the repulsive 
third virial coefficient, which sets a minimal length for collapsed blobs. 

At some value of x that corresponds to the temperature at which the microphase separa- 
tion occurs the parameter ki is always positive, so that there is no delocalization transition 
(see curve 2 in Fig. 2). Notice that assuming the existence of a maximal value for the asym- 
metry parameter, C^'^^, the reentrant localization transition is limited to Co < C™^^ > which 
imposes an inequality on x and thus on temperature. 

Notice that the parameter ki becomes infinite at some value of x- In the symmetric case 
(Co = 0) we have found that ki, ^2 become infinite at the critical value Xc = 9/(2A) (Fig.3). 
The latter means that the heteropolymer collapses onto the interface. Taking into account 
the ternary interactions will prevent this. 



III. CONCLUSION 

We have considered the adsorption of a random heteropolymer onto an interface within 
the model by Garel et al. [jll by taking into account the self-interactions between the 
monomers. The use of a preaveraging procedure within the replica method permits to map 
the present problem to a localization problem associated with a quantum mechanical Hamil- 
tonian. To find the binding state of the latter we use the variational principle based on the 
2-nd Legendre transform. We have found that self-interactions favor the localization of the 
random heteropolymer driving the delocalization transition to larger values of asymmetry. 
At the appropriate strength of self-interactions we found a reentrant localization transition 
at sufficiently high values of the asymmetry parameter Co- Although we expect that the part 
of predictions are physically irrelevant due to not taking into account the repulsive ternary 
interactions, which will effectuate a minimal collapse length, the qualitative tendencies of 
self-interactions will remain valid. The method we use can be extended to go beyond the 
approximations used in the present article. The preaveraging approximation means in fact 
a restriction of the extremum functional associated with the second Legendre transform to 
one-loop level, and can be overcomed. It is also possible to include the longitudinal behavior 
of the heteropolymer into the extremum procedure. 
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Figure captions 

Fig.l Dependence of ki (solid line) and ^2 (dashed line) as a functions of the asymmetry 
parameter at Xo = 0, Aq = 1, a = 1, T = 1 0. 

Fig. 2 Dependence of ki (solid line) and k2 (dashed line) as a functions of the asymmetry 
parameter Co at xo = 1; Aq = 1, cr = 1 (curves 1, 2 correspond to T = 0.5; 1 respectively). 
Here LS, DS correspond to localized and delocalized states. 

Fig. 3 Dependence of ki (solid line) and ^2 (dashed line) as a functions of the parameter 
X at Co = 0, Aq = 1, 0" = 1. Xc is equal to 9/ (2A). 
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Fig.l 
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